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Multiple  scattering  corrections  for  lidar  detection  of  obscured 
objects 


T.  Duracz  and  N.  J.  McCormick 


The  importance  of  multiple  scattering  for  lidar  detection  of  a  spherical  object  obscured  by  an  aerosol  is 
assessed  using  Monte  Carlo  radiative  transfer  calculations.  Multiple  scattering  correction  factors  are 
significant  and  depend  on  the  location  and  size  of  the  object,  and  the  field  of  view  and  time  resolution  of  the 
detector. 


I.  Introduction 

Lidar  detection  and  identification  of  objects  is  one  of 
the  important  tasks  in  remote  sensing.  Since  natural 
or  man-made  obscuring  aerosols  are  often  present, 
multiple  scattering  of  the  electromagnetic  radiation 
can  be  a  serious  impediment  to  effective  detection. 
Although  the  effect  of  multiple  scattering  when  char¬ 
acterizing  aerosols  has  been  studied  rather  extensive¬ 
ly,1"7  its  effect  on  lidar  returns  from  objects  is  typically 
neglected  by  assuming  a  simple  exponential  attenua¬ 
tion  along  the  optical  path. 

We  consider  the  effect  of  multiple  scattering  on 
detection  of  a  spherical  object  immersed  in  an  obscur¬ 
ing  aerosol,  as  schematically  shown  in  Fig.  1.  Both  the 
source  and  detector  are  presumed  to  be  at  the  same 
point  (so  the  lidar  is  said  to  be  monostatic).  The 
incident  lidar  pulse  at  time  t  =  0  is  located  at  (x.y.z)  = 
(0,0,0)  and  directed  along  the  x-axis.  The  space 
bounded  by  planes  x  =  0  and  x  =  x„  is  vacuum  and  the 
space  between  x  =  x„  and  x  =  x*  is  filled  with  aerosol. 
There  is  a  diffusely  reflecting  plane  boundary  at  x  =  x& 
with  albedo  at,.  For  reasons  of  symmetry,  a  diffusely 
reflecting  sphere  of  radius  r0  and  albedo  a0  immersed 
in  the  aerosol  can  be  taken  to  have  its  center  (x0ly0,z0) 
such  that 

x,  +  ra<ia<x„~r0,  y„  >  0,  z,  -  0.  (1) 

Within  the  detector  field  of  view  (FO  V)  consisting  of 
a  cone  of  half-angle  6  measured  in  radians,  all  scatter¬ 
ings  in  the  aerosol  and  reflections  from  objects  can 
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contribute  to  the  time-resolved  backscattered  signal. 
The  effect  of  multiple  scattering  increases  with  the 
FOV  half-angle  8,  which  tends  to  make  the  object  less 
discernible  even  though  the  total  received  signal  is 
greater.  The  diameter  d  of  the  detector  cone  at  the 
location  of  sphere,  for  example,  is 

d  *  2  0x„,  (2) 

and  is  greater  than  the  diameter  of  the  sphere  when 

K  >  r„.  (3) 

To  detect  and  identify  the  spherical  object,  the  fol¬ 
lowing  are  assumed  known:  distances  x0  and  xt,,  albe¬ 
do  ab,  angle  6,  and  the  single  scatter  albedo  w  and  the 
phase  function  p(f i)  of  the  aerosol,  where  /x  is  the  cosine 
of  the  angle  between  pre-  and  post-scattering  direc¬ 
tions;  since  the  radiation  is  assumed  to  remain  essen¬ 
tially  monochromatic,  the  albedo  and  phase  function 
do  not  depend  upon  wavelength.  The  phase  function 
is  normalized  such  that 

2»  j  p(n)dn  “  1.  (4) 

The  spatial  variables  x0  and  y0,  the  radius  r0,  and 
albedo  a0  of  the  spherical  object  are  the  four  un¬ 
knowns. 

The  nonanalog  Monte  Carlo  method,  which  has  been 
successfully  applied  in  other  remote  sensing  simula¬ 
tions, I>2>46  has  been  chosen  for  the  radiative  transfer 
calculations  because  of  the  configuration  of  the  3-D 
geometry.  In  this  method,  a  photon  is  tracked  from 
one  interaction  in  the  aerosol  or  with  the  object  to 
another  by  sampling  for  the  distance  to  next  interac¬ 
tion  or  escape  from  the  aerosol,  P  -  —<j~1  ln|,  where  a  is 
the  aerosol  attenuation  coefficient  in  units  of  inverse 
length  and  £  is  a  pseudorandom  number  on  the  interval 
0  <  {  £  1.  At  time  t  after  the  initial  pulse,  each  photon 
carries  a  weight  w(t )  that  is  proportional  to  the  product 
of  the  nonabsorption  probabilities  of  each  preceeding 
collision.  Whenever  there  is  an  interaction  within  the 
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reflecting 

plane 


Fig.  1.  Schematic  of  a  spherical  object  embedded  in  an  aerosol. 


detection  cone  at  time  fy,  the  it h  photon  causes  a  detec¬ 
tor  signal  Sij(t)  of 

S„(t)  -  Kwdjr- exp(-d)p(ii)i(t  -tj-  r/v),  (5) 

where  r  is  the  distance  from  the  point  of  interaction  to 
the  detector,  d  is  the  corresponding  optical  distance  to 
the  detector,  v  is  the  speed  of  light,  and  K  is  a  normaliz¬ 
ing  factor  that  is  independent  of  time.  The  contribu¬ 
tions  to  the  signal  from  all  scattering  and  reflections 
from  the  object  or  background  surface  are  summed  and 
placed  in  appropriate  time  bins  to  obtain  the  time- 
resolved  simulated  backscattered  signal. 

Before  analyzing  the  effects  of  multiple  scattering  it 
is  instructive  to  consider  the  case  of  an  unobscured 
spherical  object  (or  equivalently,  one  obscured  by  a 
nonabsorbing  aerosol  that  is  a  perfect  straight-ahead 
scatterer).  If  an  instantaneous  pulse  is  emitted  at 
time  t-  0  in  the  direction  of  the  * -axis  and  the  spheri¬ 
cal  object  is  centered  at  (zatyo,0),  the  total  returned 
signal  received  by  the  detector  is  given  by 

S(C)  -  K(aJ2r)r~2Ht  -  2r/u),  y0  <  r0. 


either  the  aerosol  or  the  object.  If  the  aerosol  is  opti¬ 
cally  thin  enough  or  the  detector  FOV  angle  d  is  small 
enough,  only  once-scattered  photons  need  be  consid¬ 
ered.  The  returned  signal  Stl)(t)  received  by  the  de¬ 
tector  is  then 

S"’(f)  «  K  exp[-2d(t/2)]l2(vtr-wp(-l) 

X  feu  -  2x„/v)  -  eu  -  2r/u)] 

+  r~'(aJ2x)i(t  -  2r/o)|,  y0  <  r0, 

“  K  eip[-2d(t/2)]2(ut)'i«p(-l)e(t  -  2xju),  y„  >  r„,  (10) 

where  d(t)  =  a(vt  -  xa )  is  the  optical  distance  from  an 
interaction  point  to  the  detector,  r  is  the  distance  from 
the  object  to  the  detector,  and  9(t)  is  the  step  function: 

e(t)  =  1,  t  >  o, 

=  0,  t<0.  (11) 

When  the  aerosol  is  so  thick  that  multiple  scatter¬ 
ings  must  be  taken  into  account,  the  toted  signal  re¬ 
ceived  by  the  detector  from  all  the  photons  is 

S(t)  -  £  =  B(t)  +  R{t)’  U2) 

which  is  the  sum  of  the  backscattered  signal  B(t)  gen¬ 
erated  by  one  or  more  scatterings  in  the  aerosol  and  the 
reflected  signal  f?(t)  from  at  least  one  reflection  from 
the  object,  with  possibly  one  or  more  scatterings  before 
or  after  reflection.  In  a  sense,  B(t )  represents  pure 
noise  since  there  is  no  information  about  the  object. 

The  signal  obtained  in  an  actual  measurement  is 
always  averaged  over  a  small  interval  of  time.  Like¬ 
wise,  in  the  Monte  Carlo  method  this  occurs  when  the 
detector  count  rate  is  tallied  in  the  fcth  time  bin  A t*. 
In  practice  this  means  that  S(t)  in  Eq.  (12)  is  of  the 
form 


» 0,  otherwise.  (6) 

This  result  follows  from  Eq.  (5)  with  d  =  0  and  w  *  1, 
since  there  is  only  a  single  reflection  from  the  object 
and  no  attenuation  along  the  path,  and  since  p(n)  =  a  J 
2r  corresponding  to  an  equal  probability  for  scattering 
in  any  direction  leaving  the  surface.  From  the  time 
delay  td  =  2 r/o  after  the  initial  pulse  and  the  magni¬ 
tude  of  the  time-integrated  return  S0  =  Kr~2a0/2r  in 
Eq.  (6)  it  is  possible  to  estimate  the  distance  to  the 
sphereandits“albedo7 

r-  tdul 2.  (7) 

a„  -  SQx(tdu)2/2K.  (8) 

The  location  (x0y0)  of  the  object  and  its  radius  r0  are 
linked  to  the  distance  r  through 

<x,-r)2  +  ^-rj.  (9) 

It  is  worth  noting  that  the  localized  return  in  time  in 
Eq.  (6)  depends  on  a  localized  spot  illumination.  With 
a  larger  area  of  illumination,  however,  there  can  be 
considerable  geometric  pulse  broadening.8 

When  an  aerosol  is  present,  Eq.  (6)  must  be  modified 
to  account  for  the  possibility  of  backscattering  from 


Sv  -  v  f  SJt)dt,  (13) 

where  the  time  interval  At*  is  defined  such  that 

ar*  -  [(A  -  i)At,  *At],  k  >  i.  (14) 

Multiple  scattering  can  be  expressed  in  the  form  of  a 
multiple  scattering  correction, 


which  is  the  ratio  of  the  total  return  to  the  single 
scattering  return  from  Eq.  (10),  with  N  the  number  of 
photons  used  in  the  simulation.  The  multiply-scat¬ 
tered  returned  signal  S(t)  —  SU)(t)  is  noncoherent  due 
to  the  random  nature  of  the  scatterings  in  the  aerosol. 
Thus  this  portion  of  the  returned  signal  would  contrib¬ 
ute  directly  to  the  total  signal  measured  by  a  noncoher¬ 
ent  lidar  detector,  as  assumed  in  this  paper,  or  would 
contribute  to  the  background  noise  of  a  coherent  (het¬ 
erodyne)  lidar  detector. 

fl.  Procedure  for  Numerical  Tests 

The  phase  function  pg(m)  of  the  obscuring  aerosol 
was  obtained  from  a  Mie  scattering  calculation  for 
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spherical  particulates  with  size  parameter  a  =  2ira/\, 
where  a  is  the  sphere  radius  and  X  is  the  radiation 
wavelength.  A  Gaussian  distribution  of  sizes  was  as¬ 
sumed8  with  a  mean  value  of  (o)  =  10  and  a  standard 
deviation  of  <r  =  2.  A  nonabsorbing  medium  (to  =  1) 
with  a  refractive  index  of  n  -  ik  =  1.3  was  selected,  and 
this  resulted  in  a  moderately  forward  oriented  scatter¬ 
ing.8 

The  unit  of  length  for  the  calculations  was  taken  to 
be  the  mean  free  path  between  scatterings  in  the  aero¬ 
sol,  and  the  unit  of  time  was  the  mean  time  between 
scatterings,  so  that  v  =  1.  The  distance  from  the  lidar 
to  the  edge  of  the  aerosol  was  arbitrarily  set  to  xa  =  3; 


Fig.  2.  Returned  signal  vs  time  for  a  sphere  of  radius  0.01  with  its 
center  at  distance  4  from  a  lidar  whose  FOV  half-angle  is  5  X  10-3 
rad.  Time  bin  duration  is  0.05. 


•»i 


Fig.  3.  Returned  signal  normalized  to  that  in  Fig.  2  vs  time  for  a 
sphere  of  radius  0.01  with  its  center  at  distance  5  from  a  lidar  whose 
FOV  half-angle  is  5  X  10~3  rad.  Time  bin  duration  is  0.05. 


this  distance  translates  the  time  axis  of  all  calculations 
and  also  influences  the  geometry  of  the  problem. 

An  important  parameter  in  the  study  was  the  time 
resolution  of  the  detector,  i.e.,  the  duration  of  each 
time  bin  At.  The  values  varied  from  0.05  to  0.25. 

To  maximize  the  lidar  return,  a  nonabsorbing  sphere 
(a0  =  1)  was  placed  on  the  x-axis  (y0  -  0).  Two  values 
of  the  sphere  radius  were  selected,  r0  =  0.01  and  0.1, 
while  the  background  plane  boundary  was  placed  far 
away  from  the  sphere  (xt>  =  « ).  The  values  of  x0  and  6 
remained  variable,  in  addition  to  At,  so  S*  = 
Sk(x0,9,At)  and  C*  =  C*(xo,0,At). 

An  estimate  of  the  relative  error  e*  associated  with 


Fig.  4.  Returned  signal  normalized  to  that  in  Fig.  2  vs  time  for  a 
sphere  of  radius  0.01  with  its  center  at  distance  5  from  a  lidar  whose 
FOV  half-angle  is  5  X  10“-  rad.  Time  bin  duration  is  0.05. 


Fig.  5.  Returned  signal  vs  time  for  t  sphere  of  radius  0.01  with  its 
center  at  distance  4  from  a  lidar  whose  FOV  half-angle  is  5  x  10~3 
rad.  Time  bin  duration  is  0.25. 
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Fig.  6.  Returned  signal  normalized  to  that  in  Fig.  5  vs  time  for  a 
sphere  of  radius  0.01  with  its  center  at  distance  5  from  a  lidar  whose 
FOV  half-angle  is  5  X  10"3  rad.  Time  bin  duration  is  0.25. 


Si,  for  N  photon  histories  can  be  obtained  from 


For  large  N  this  becomes  identical  with  the  estimated 
standard  deviation  of  the  mean  tallied  signal  "divided 
by  the  mean  signal.  The  computer  simulations  were 
made  with  N  -  5  X  105  or  106,  and  the  value  of  e*  did 
not  exceed  a  few  percent  even  for  cases  with  the  small¬ 
est  time  bins  At  =  0.05. 

<11.  Effect  of  Multiple  Scattering 

Examples  of  the  simulated  signal  return  S*  obtained 
as  in  Eq.  (13)  are  shown  in  Figs.  2-4  for  a  sphere  of 
radius  r0  -  0.01  and  time  bins  of  duration  A  =  0.05. 
The  values  are  normalized  to  the  maximum  value  of 
unity  in  Fig.  2. 

In  Fig.  2  the  distance  to  the  lidar  is  x„  -  4  and  the 
detector  FOV  half-angle  is  6  -  5  X  10-3  rad,  while  in 
Fig.  3  the  distance  to  the  lidar  has  been  increased  to  x0 
■  5  with  the  same  FOV.  The  strong  damping  effect  of 
the  additional  intervening  aerosol  is  visible.  "lb  Fig.  4 
the  distance  is  x0  =  5,  as  in  Fig.  3,  but  the  FOV  has  been 
increased  to  8  =  5  X  10"2  rad.  The  received  signal  is 
stronger  because  of  additional  radiation  scattered 
from  the  object  and  aerosol  and  a  slight  spreading  in 
time  of  the  signal  also  occurs. 

In  Figs.  5-7  all  the  parameters  are  identical  to  those 
in  Figs.  2—4  except  that  At  m  0.25  instead  of  0.05.  As 
At  increases  the  peak  value  of  the  returned  signal 
appears  smaller  relative  to  the  returned  signal  in  other 
bins  because  of  the  aerosol.  Also,  by  comparing  Figs.  4 
and  7  an  apparent  broadening  of  the  returned  signal 
near  the  peak  can  be  seen  for  the  case  with  larger  At. 
This  broadening  is  sensitive  to  the  precise  location  of 
the  At/,  relative  to  twice  the  transit  time  between  the 


Fig.  7.  Returned  signal  normalized  to  that  in  Fig.  5  vs  time  for  a 
sphere  of  radius  0.01  with  its  center  at  distance  5  from  a  lidar  whose 
FOV  half-angle  is  5  X  10~2  rad.  Time  bin  duration  is  0.25. 


lidar  and  the  object,  2(x0  -  r0),  since  it  may  happen 
that  part  of  the  contribution  of  the  returned  signal 
from  the  object  gets  recorded  in  adjacent  time  bins. 

Because  the  peak  value  of  the  returned  signal  is  of 
principal  interest,  we  now  focus  on  the  correction  fac¬ 
tor  Csp(x,e,At),  the  Ck(x0,8,At)  in  Eq.  (15)  for  a  time  bin 
A t,p  symmetrically  located  about  the  time  of  the  peak 
signal  return: 

At,„  -  [2(x„  -  r„)  -  At/2.  2(x„  -  r„)  +  At/2].  (17) 

C,p  -  1  in  the  limiting  case  of  no  aerosol  (x0  =  x0)  and 
when  multiple  scattering  is  insignificant  (0  -*  0  or  At 
— •  0).  For  the  values  of  our  example  calculations,  the 
principal  changes  in  Cjp  occur  for  x,  <  x,  i  x,  +  2 
(where  x0  =  3)  and  10-3  <  9  <  10_1  rad. 

The  results  for  C,p  with  a  spherical  object  of  radius  r0 
-  0.01  and  At  =  0.05, 0.15,  and  0.25  are  shown  in  Figs. 
8-10.  (Note  the  different  scales  for  the  Csp  values.) 
The  multiple  scattering  correction  increases  to  signifi- 


Fig.  S.  Multiple  scattering  correction  factor  for  the  symmetric  peak 
time  bin  for  a  sphere  of  radius  0.01  vs  the  lidar-to-sphere  center 
distance  and  detector  FOV  half-angle.  Symmetric  time  bin  dura¬ 
tion  ia  0.05. 
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5.0  -3.0 

Fig.  9.  Multiple  scattering  correction  factor  for  the  symmetric  peak 
time  bin  for  a  sphere  of  radius  0.01  vs  the  lidar-to-sphere  center 
distance  and  detector  FOV  half-angle.  Symmetric  time  bin  dura¬ 
tion  is  0.15. 


5.0  -3.0 

Fig.  10.  Multiple  scattering  correction  factor  for  the  symmetric 
peak  time  bin  for  a  sphere  of  radius  0.01  vs  the  lidar-to-sphere  center 
distance  and  detector  FOV  half-angle.  Symmetric  time  bin  dura¬ 
tion  is  0.25. 


5.0  -3.0 

Fig.  11.  Multiple  scattering  correction  factor  for  the  symmetric 
peak  time  bin  for  a  sphere  of  radius  0.1  vs  the  lidar-to-sphere  center 
distance  and  detector  FOV  half-angle.  Symmetric  time  bin  dura¬ 
tion  is  0.05. 


cant  values  for  larger  distances  x„  and  larger  FOV  half- 
angles  6.  The  effect  is  even  more  pronounced  for  a 
sphere  of  radius  r0  ■  0.1,  as  seen  in  Figs.  11  and  12  for 
At  -  0.05  and  0.25,  respectively. 

The  rising  values  of  C,p  also  reflect  the  fact  that  the 
returned  signal  decreases  rapidly  as  the  lidar-to- 
sphere  distance  increases.  Figure  13  shows  the  de¬ 
crease  in  S,p,  the  normalized  returned  signal  in  bin 
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Fig.  12.  Multiple  scattering  correction  factor  for  the  symmetric 
peak  time  bin  for  a  sphere  of  radius  0.1  vs  the  lidar-to-sphere  center 
distance  and  detector  FOV  half-angle.  Symmetric  time  bin  dura¬ 
tion  is  0.25. 


5.0  -3.0 

Fig.  13.  Normalized  returned  signal  for  a  sphere  of  radius  0.01  vs 
the  iidar-to-sphere  center  distance  and  detector  FOV  half-angle. 
Symmetric  time  bin  duration  is  0.25. 


Fig.  14.  Multiple  scattering  correction  factor  for  a  fixed  time  bin 
grid  for  a  sphere  of  radius  0.01  vs  the  lidar-to-sphere  center  distance 
and  detector  FOV  half-angle.  Time  bin  duration  is  0.25. 


At ,p,  for  r0  *  0.01  and  At  *  0.25.  A  very  similar  pattern 
occurs  for  any  other  of  the  assumed  r0  and  At  values. 

If  the  multiple  scattering  correction  factor  is  calcu¬ 
lated  using  a  fixed  time  bin  grid,  as  in  Eq.  (14),  the  peak 
return  signal  5(f)  will  sometimes  occur  near  the  begin¬ 
ning  or  end  of  a  time  bin,  as  mentioned  earlier,  and  a 
fraction  of  the  time-broadened  peak  return  will  be 


Fig.  15.  Multiple  scattering  correction  factor  for  a  fixed  time  bin 
grid  for  a  sphere  of  radius  0.01  vs  the  lidar-to-sphere  center  distance. 
The  detector  FOV  half-angle  is  10"1  and  the  time  bin  duration  is  0.25 
rad. 

counted  in  an  adjacent  bin.  Figure  14  illustrates  the 
correction  factor  C/p  for  such  a  case  of  a  /ixed  time  bin 
grid  with  At  -  0.25  and  r0  =  0.01;  Cfp  is  the  same  as  C»  of 
Eq.  (15)  for  the  time  bin  k  in  which  the  peak  occurs. 
The  figure  should  be  compared  with  Fig.  10  for  Csp  that 
was  calculated  for  the  symmetric  time  bin  in  Eq.  (17). 
There  is  a  clear  periodic  pattern  in  Fig.  14  that  appears 
even  more  dramatic  in  the  enhanced  view  in  Fig.  15  for 
the  case  of  the  large  FOV  half-angle  of  0  -  10~l  rad. 
The  spacing  of  the  pattern  corresponds  to  2A t,  the 
additional  transit  time  for  the  signal  if  the  sphere  is 
near  the  end  of  the  time  bin  in  the  fixed  grid  instead  of 
at  the  beginning. 

IV.  Conclusions 

Monte  Carlo  calculations  have  been  performed  to 
assess  the  effect  of  multiple  scattering  on  a  lidar  return 


signal  from  a  spherical  object  obscured  by  an  aerosol. 
Figures  2-7  demonstrate  the  expected  sharp  peak  in 
the  returned  signal  caused  by  the  presence  of  the  ob¬ 
ject,  with  some  signal  decrease  and  broadening  for 
larger  distances  and  detector  FOVs.  The  importance 
of  multiple  scattering  is  illustrated  in  Figs.  8-12  for  a 
time  bin  symmetrically  located  about  the  peak  value  of 
the  time-dependent  return  signal.  W e  found  that  it  is 
important  to  distinguish  between  the  use  of  a  symmet¬ 
ric  time  bin  located  about  the  peak  signal  return  (as  in 
Fig.  10)  and  a  fixed  bin  location  (as  in  Fig.  14)  in 
calculating  the  multiple  scattering  correction  factor. 
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